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^ ' We calculate the one-loop non-factorizable QCD corrections to the production and decay of 

pairs of top quarks at various collider experiments. These non-factorizable corrections inter- 
connect the different production and decay stages of the off-shell top-pair production processes. 
This in particular affects the invariant-mass distributions of the off-shell top quarks, resulting 



P i' 

Q I in a shift of the maximum of the distorted Breit-Wigner distributions. Although the non- 

i-^ ' factorizable corrections can be large, the actual shift in the mass as determined from the peak 

position of the corrected Breit-Wigner line-shape is below 100 MeV. 
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1 Introduction 

At present and future collider experiments, a detailed investigation of the production of top- 
quark pairs will substantially contribute to our knowledge of the top-quark properties and 
thereby of the Standard Model. An improved measurement of the top-quark mass rrit, for 
instance, can serve to obtain improved indirect sensitivity to the mass of the Standard Model 
Higgs boson. This is achieved by combining the high-precision measurements of the electroweak 
parameters at LEP/SLC with the direct measurements of the top-quark and IV-boson masses. 

Pairs of top quarks can be produced in hadron collisions at the Tevatron {pp) and LHC {pp), 
as well as in e~^e~ and 77 collisions at a future linear collider. Since the top quark has a large 
width as compared to the QCD hadronization scale, Tt ~ 1.4 GeV » Aqcd ~ 200 — 300 MeV, 
it predominantly decays before hadronization takes place. Therefore the perturbative approach 
can be used for describing top quarks. The main lowest-order (partonic) mechanisms for the 
pair production of top quarks are 

e"'"e~, 77 — > tt -^ hW^hW^ -^ 6 fermions, (1) 

QQi 99 ^ tt ~^ bW^bW — > 6 fermions. (2) 

A lot of effort has been put into an adequate theoretical description of these reactions (see 
e.g. Ref. for two review papers). Most of these studies treat the top quarks as stable particles, 
which is a reasonable approximation since Tt/rrit = (9(1%). For the reactions qq,99 — > tt these 
studies comprise QCD [§] and electroweak [Q one-loop corrections, as well as the resummation of 
soft-gluon effects |^. Also for the reactions e"'"e~, 77 -^ tt both the QCD [Q and electroweak 
one-loop corrections are known. Moreover, the ti threshold, with its sizeable QCD |j^ and 
Yukawa interactions ||^, has been analysed in detail. 

One would, however, like to treat the top quark as an unstable particle, with a Breit-Wigner 
distribution describing its line shape. The most economic approach for treating processes that 
involve the production and subsequent decay of unstable particles is the so-called leading- 
pole approximation (LPA) 0]. This approximation is based on an expansion of the complete 
amplitude around the poles of the unstable particles, which can be viewed as a prescription for 
performing an effective expansion in powers of Fj/Mj. Here Mj and Fj stand for the masses 
and widths of the various unstable particles. The residues in the pole expansion are physically 
observable and therefore gauge- invariant. The actual approximation consists in retaining only 
the terms with the highest degree of resonance. In the case of top-quark pair production only the 
double-pole residues are hence considered and the LPA becomes a double-pole approximation 
(DPA). This approximation will be valid sufficiently far above the tt-threshold. If in reactions 
(|I|) and (0) also the W bosons are treated as unstable particles, then also for these particles 
the leading pole residues should be taken. In this approach the complete set of corrections 
to reactions (^ and (H) naturally sphts into two groups: factorizable and non-factorizable 
corrections. The factorizable corrections are directly linked to the density matrices for on- 
shell production and decay of the unstable particles. The non-factorizable corrections can be 
viewed as describing interactions that interconnect different (production/decay) stages of the 
off-shell process. A detailed discussion of this method with all its subtleties can be found in 



Ref. [10|, where the method has been apphed to the complete set of 0{a) radiative corrections 



to the process e~^e — > W^W -^ 4 fermions. For tt production partial results along this 



line exist |T^, involving a subset of the factorizable corrections to the reaction e+e~ — > tt ^ 
bW^bW~ — > 6 fermions. However, the non-factorizable corrections are needed for a complete 
0{as) calculation. 

In recent years the necessary methods for calculating such non-factorizable corrections have 
been developed. In Ref. [1^1 a first complete calculation was performed for tt production in e~^e~ 



collisions. In Ref. [|I3] the non-factorizable corrections were calculated for W-pair production, 
revealing differences with the results of Ref. |jl2| . The results of Ref. |jl3l were confirmed by an 



independent calculation |T^ as well as by a re-analysis of the results of Ref. |]T2 



In this paper we apply our calculations presented in Ref. |13| to the non-factorizable 0{as) 
corrections to ti production at various colliders. We discuss the effect on the invariant-mass 
distribution of the off-shell top quark and the resulting shift in the maximum of the distorted 
Breit-Wigner distribution. 



2 Definition of the non-factorizable corrections 

In the LPA approach reactions like (|l]) and (^, which involve unstable particles during inter- 
mediate stages, can be viewed as consisting of separate subprocesses, i.e. the production and 
decay of the unstable particles. Having this picture in mind, the complete set of radiative 
corrections can be separated naturally into a sum of corrections to these subprocesses, called 
factorizable corrections, and those corrections that interconnect various subprocesses, called 
non-factorizable corrections. It should be noted, however, that it is often misleading to identify 
the non-factorizable contributions on the basis of diagrams. Such a definition is in general 
not gauge- invariant. Rather one should realize that only real/virtual semi-soft gluonsQ with 
Eg = 0(Ti) will contribute, the contributions of the hard gluons being suppressed by Vi/ Eg. 
This is a consequence of the fact that the various subprocesses are typically separated by a big 
space-time interval of Oil/Vi) due to the propagation of the unstable particles. The subpro- 
cesses can be interconnected only by the radiation of semi-soft gluons with energy of 0{Ti), 
which induce interactions that are sufficiently long range. Hard gluons [Eg = 0{Mi) ^ Fj) as 
well as massive particles induce short-range interactions and therefore contribute exclusively to 
the factorizable corrections, which are governed by the relatively short time interval ~ 1/Mj on 
which the decay and production subprocesses occur. A more detailed discussion of these issues 
can be found in Refs. [llOl |T5|, HI. 



In Fig. [^ we show schematically the partonic process qq ^ tt ^>- bW~^bW~ —>■ 6 fermions. 
The process consists of five subprocesses, which we will denote by ttprod, ^dec, ^dec, W^dec ^^^ 
W^^. In Fig. m these subprocesses are indicated by the open circles. The non-factorizable 
semi-soft gluon interactions interconnect any two different subprocesses, as is exemplified in 

^These gluons will still be perturbative in our case as their typical energy {Eg ~ Tt^w ^ 1.4 GeV) largely 
exceeds the QCD hadronization scale (Aqcjj w 200 — 300 MeV). 




Figure 1: The generic structure of the complete tt-production process qq ^ ti ^ bW~^bW~ — > 
6 fermions in the LPA. The open circles denote the various production and decay subprocesses. 
As an example also the non-factorizable semi-soft gluon interaction between the two top-quark 
decay subprocesses is shown. 

Fig. |T| for the two top-quark decay subprocesses. The coupling of such a gluon to a certain 
subprocess can be written in terms of semi-soft currents. In contrast to soft-gluon currents, the 
effect of the gluon momentum on the unstable-particle propagators cannot be neglected in the 
semi-soft currents. The various non-factorizable corrections to the cross-section are just given 
by all possible interferences of the semi-soft currents. This will be made more explicit in the 
next section. 



3 Colour dependence of the non-factorizable corrections 

We start off by considering the simpler case of stable W bosons. At the end of this section 
we will indicate what happens if the W bosons decay hadronically. For stable W bosons one 
can identify three subprocesses: ttprod, ^dec, ^dec- The non-factorizable corrections are given 
by the semi-soft gluon interferences between these different subprocesses. As only semi-soft 
gluons contribute, the virtual and real matrix elements factorize in terms of lowest-order matrix 
elements and semi-soft currents. In view of the possible presence of coloured particles in the 
initial state {qq,gg), this factorization depends on the colour structure. For the reactions (|l]), 
which involve only colourless initial-state particles, the ti pair is produced in a singlet state. 
In contrast, the ti pair is produced in an octet state in the lowest-order annihilation process 
qq —>■ ti, which involves the time-like exchange of a gluon. Both singlet and octet states are 
present in the lowest-order gluon-fusion reaction gg — > ti, since in that case also space-like 
top-quark-exchange diagrams contribute. Because of these differences in the colour structure 
of the lowest-order reactions, also the non-factorizable corrections will come out differently, as 
we will see from the following discussion. 

In order to keep the notation as general as possible, we write the lowest-order partonic 
reactions in the generic form 

Qiiqi)Q2iq2) ^ tip^)iip,) ^ bih)W+ik[) bih)W-ik',), (3) 



where Q1Q2 = {e'''e~, 77, 55,(7(7}. The corresponding lowest-order matrix element will be de- 
noted by (A^o)ij'^S where i,j indicate the t, i colour indices in the fundamental representation. 
The colour indices Ci, C2 belonging to Qi, Q2 depend on the specific initial state: they are absent 
for the colourless e^e~ and 77 initial states, and they are in the fundamental/adjoint repre- 
sentation for the qq/ gg initial states. The momentum, Lorentz index, and colour index of the 
semi-soft gluon will be denoted by A;, /i, and a, respectively. 

By using the relation 
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Ikl 
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(4) 



for the SU{N) generators T" in the fundamental representation (with N = 3 for QCD), the 
virtual and real non-factorizable corrections take the generic form: 
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Here the pre-factor consists of the lowest-order phase-space factor in the LPA [dVf^ , the partonic 
flux factor [l/(2s)], and the initial-state spin and colour average [l/i^jn]. The non-factorizable 
kernels can be expressed in terms of semi-soft currents according to 
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The terms proportional to diUjudjiii project on the lowest-order singlet tt states, whereas the 
terms proportional to dimi Sj/j" completely factorize the lowest-order cross-section. The colour 



structure 



for e'^e ,77 



Qtn = { '^«(^")4'4 + '^4'4(^")c;<' for qq 
'^<c;'(^")<4 +4^'4(^")c;c'/ for gg 



(9) 



depends on the specific initial state and in general does not project on explicit lowest-order 
ti colour states. Here F" are the SU{N) generators in the adjoint representation, which are 
defined in terms of the SU{N) structure constant according to {F"-)bc = —if"'^'^. Note that for 
e"'"e~ and 77 initial states the currents jTg, Jj^ and Xq completely drop out of Eqs. (^ and (^, 
as it should be for colourless particles in the initial state. 

The semi-soft currents appearing in the virtual non-factorizable corrections are given by 
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for gluon emission from the decay stages of the process, and 
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(11) 



for gluon emission from the production stage of the process. Here gs is the QCD gauge coupling 
and D12 = P12 ~ ^ + i^^t^t is a shorthand notation for the inverse top-quark propagators. 
Note the difference in the sign of the io parts appearing in the currents J@,J@ and Jq.Jq. 
These infinitesimal imaginary parts are needed to ensure a proper incorporation of causality. 



The corresponding semi-soft real-gluon currents read 
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By simple power counting one can explicitly see from the above specified currents that 
the contributions of hard gluons are suppressed and that effectively only semi-soft gluons with 
Eg = ko = 0(Tt) contribute. In view of the pole structure of the virtual corrections, governed by 



the infinitesimal imaginary parts io, many of the non-factorizable corrections will vanish when 
virtual and real-gluon corrections are added up. For instance, all initial-final state interferences 
will vanish, leaving behind a very limited subset of 'final-state' interferences [jlOl , |T6|. The 
following holds for the remaining interferences: 
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with similar effective replacements for Jii. As a result of these properties of the non-factorizable 
corrections, a factorization per colour structure emerges: 
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Here (icTBom,! and (icrBorn,8 are the lowest-order multi-differential cross-sections for producing the 
intermediate ti pair in a singlet and octet state, respectively. For completeness we note that 
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(16) 



The non-factorizable factor 6n{ can be obtained from Ref. ||13[. The results of Section 4 of that 
paper should be used, since those allow for massive decay products from the unstable particles, 
which is the case for the top-quark decay. 

We conclude by considering the case that also the W bosons are unstable. This adds two 
decay subprocesses, VF^c ^^^ ^dec ^^ ^^^ three we have considered so far. If the W bosons 
decay leptonically, nothing changes as the gluon cannot couple to the W decay subprocesses 
in that case. For a hadronically decaying W boson additional interferences have to be taken 
into account. However, such interferences trivially vanish as a result of the singlet nature of 
the VT-boson decays [i.e. Tr(T")=0]. 



4 Numerical results 



With the help of Eqs. (|I^)-([T^) we can now in principle evaluate all kinds of multi-differential 
distributions, with and without non-factorizable corrections. Although the factorized structure 
of the non-factorizable corrections is very transparent in Eq. (pi^D, integration of the multi- 
differential cross-sections will affect this structure. For instance, in Eq. (|1^ the correction to 
the singlet cross-section differs by a factor —8 with respect to the octet one. However, for 
the calculation of the relative non-factorizable corrections to a one- dimensional distribution, 
one has to evaluate the ratio of the integrated Eqs. ([l^) and ([T6|) . Since 6n{ depends on the 
integration variables, the thus-obtained singlet and octet correction factors will not necessarily 
differ by the factor —8. 



6 



At this point we stress that any observable that is inclusive in both top-quark invariant 
masses, such as the total cross-section, will not receive any non-factorizable corrections. This 



is a typical feature of these interconnection effects [|I^ . As an example of a distribution that is 
subject to non-vanishing non-factorizable corrections we focus on the invariant-mass distribu- 
tion of the top quark, which can be used for the mass determination. To this end we determine 
the non-factorizable correction 6n{{M) for the distribution 

da daBorn 



dM dM 



1 + 5„f(M) , (17) 



where M is the invariant mass of the 6-quark and the W~^ boson. The maximum of the Breit- 
Wigner distribution can be used to determine the top-quark mass. The linearized shift of this 
maximum as induced by the non-factorizable corrections is given by 



AM = - r^ ^ ' 



dM 



M=mt 



The correction 6^f{M) is calculated for the four different mechanisms of ti production, 
i.e. initiated by e~^e~, 77, qq and gg. For the centre-of-mass energies of these (partonic) 
reactions we take ^/s = 355 GeV and 500 GeV. These values exemplify the non-factorizable 
corrections in the vicinity of the threshold and far above it. As mentioned before, the adopted 
approximation in our calculation (LPA) forces us to stay sufficiently far above the ti threshold 
(read: a few times r^). The numerical values for the input parameters are 

mt = 173.8 GeV, Mw = 80.26 GeV, Mz = 91.187 GeV, (19) 

and 

ri = 1.3901 GeV, (20) 

the latter being the 0{as) corrected top-quark width. The correction S^f is proportional to 
as, for which we have to choose the relevant scale. For y/s = 355 GeV the main contribution 
originates from the non-factorizable Coulomb effect present in 5nf. Its typical momentum is 
determined by the top-quark width Tt and velocity /?: Tt/P ~ 6.8 GeV. At 500 GeV softer 
gluons contribute and therefore the typical gluon momentum is F^ ~ 1.4 GeV. Therefore we 
choose 

as(l-4 GeV) ^ 0.3536 for v^ = 500 GeV, (21) 

as(6.8 GeV) ^ 0.1955 for v^ = 355 GeV, (22) 

corresponding to as{Mz) = 0.1180 at the Z peak. It should be noted that choosing another 
scale in as will only affect the normalization of the correction. 

In Fig. ^ the non-factorizable correction (5nf is plotted as a function of the invariant mass 
M at the centre-of-mass energy of 355 GeV. The ^nf values for the pure singlet e^e~ initial 
state and the pure octet qq initial state differ approximately by the afore-mentioned factor of 
—8. For the gg initial state the Born octet part is larger than the singlet one, resulting in a 
non-factorizable correction that lies between the e+e~ and the qq case. The correction for the 
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Figure 2: The relative non-factorizable correction 5„j(M) to the single invariant-mass distri- 
bution da/dM. Centre-of-mass energy: y/s = 355 GeV. 
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Figure 3: The relative non-factorizable correction Snf{M) to the single invariant-mass distri- 
bution da/dM. Centre-of-mass energy: ^i = 500 GeV. 



77 initial state is virtually indistinguishable from the e^e'^ one and is therefore not displayed. 
Evidently the distortion effects from the singlet corrections are very large, which is due to a large 
non-factorizable Coulomb correction inside 5nf. The maximum of the Breit-Wigner distribution 
is hardly affected by this large correction. One finds for the various initial states e^e~(77), gg 
and qq AM ^ —85, —15 and -l-lOMeV respectively. The situation at 500 GeV is depicted in 
Fig. 1^. The overall correction is small, which is typical for non-factorizable corrections further 
away from threshold. The shift in the maximum of the Breit-Wigner distribution is of the 
order of 5 MeV for the e^e~ and 77 initial states, and even smaller for the qq and gg initial 
states. 

In order to obtain hadronic distributions from the partonic ones, the results for the qq and 
gg initial states should of course be properly folded with the parton densities of the colliding 
hadrons {pp at the Tevatron, pp at the LHC). The bulk of the partonic contributions originates 
from the energy region not far above the tt threshold (s < Sm^ , i.e. y/s < 500 GeV), which is 
exemplified by the partonic energies 355 and 500 GeV used in our analysis. 



5 Conclusions 

In this paper we have summarized the gauge-invariant description for calculating the 0{as) 
non-factorizable QCD corrections to pair production of top quarks. The formalism is presented 
in a general way, making it applicable to all relevant initial states. The resulting final formula for 
the non-factorizable corrections involves the same quantity S^t for all reactions. This quantity 
can be numerically calculated using expressions available in the literature. 

Although the formalism can be used for numerical studies of many distributions, the focus 
of our numerical evaluation has been on the invariant-mass distribution of the top quark, which 
can be used for extracting the top-quark mass. In spite of the possible sizeable deformations 
of this line-shape distribution, its maximum is shifted by less then 100 MeV. Therefore, if 
the top-quark mass is extracted experimentally from the peak position of the line-shape, the 
non-factorizable corrections can be safely neglected. If the precise shape of the Breit-Wigner 
distribution is used in the experimental analysis, the non-factorizable corrections should be 
taken into account properly. In particular if the singlet colour state dominates. In addition 
higher-order non-factorizable corrections might be needed. 
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